NASA TN d-1562 


N 


NASA TN D-1562 



TECHNICAL NOTE 

-Crp D-1562 . j 


METHOD FOR DESIGN OF PUMP IMPELLERS USING A 

HIGH-SPEED DIGITAL COMPUTER 

By Norbert O. Stockman and John L. Kramer 

Lewis Research Center 
Cleveland, Ohio 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON January 1963 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


TECHNICAL NOTE D-1562 


METHOD FOR DESIGN OE PUMP IMPELLERS USING A 
HIGH-SPEED DIGITAL COMPUTER 
By Norbert 0. Stockman and John L. Kramer 


SUMMART 

A method of designing pump impellers is derived from the equations of motion 
and continuity for incompressible nonviscous relative flow. The flow is assumed 
to follow a known stream surface (representing blade shape) that extends from hub 
to shroud. Equations are also derived for approximate blade-surface velocities 
and pressures. A detailed numerical procedure and block diagram are given for 
use on a digital computer. A numerical example that illustrates limited use of 
the method is presented and further uses are indicated. 


INTRODUCTION 

Centrifugal pumps are being considered for space applications in chemical 
and nuclear-rocket engines using cryogenic fluids and in power- conversion systems 
using liquid metals. For these applications, the pumps should be of minimum size 
and weight, should have high efficiency, and should operate cavitation-free or 
with controlled cavitation at high rotative speeds. A knowledge and control of 
the internal flow are necessary to meet these requirements. 

A method based on stream- filament techniques for analyzing the flow in the 
impeller of a centrifugal air compressor was developed at the Lewis Research 
Center (ref. 1). This method determines velocities, pressures, and streamlines 
in the meridional plane (a plane containing the impeller axis) and of approximate 
velocities and pressures on the blade surfaces. From the analysis an impeller 
design method (ref. 2) was derived that allowed direct control of the hub veloc- 
ities and indirect control of the velocities and pressures in the impeller. This 
method enabled the designer to avoid velocity gradients conducive to boundary- 
layer separation and to avoid eddy formation on the blade surfaces. 

Tie method of reference 2 was used to redesign the shroud shapes of three 
existing impellers for air compressors for which experimental data were avail- 
able. Data for the redesigned impellers were obtained, and the comparison of 
these two sets of data showed significant gains in performance of all three 
impellers ( ref. 3 ) . 


These favorable results with air compressors, together with recent pump 
applications that require knowledge and control of the internal flow, led to the 


use of this method for the design of centrifugal pumps. To facilitate this use, 
it is desirable to determine the specific flow equations and design methods for 
the incompressible flow in pump impellers, and to detail the procedures for high- 
speed digital- computer use. This report presents the development of such a 
method. 

The basic flow equations are derived from the equations of motion for a 
mean-flow surface that extends from impeller hub to shroud. The required equa- 
tion for stream tube volume flow is derived from continuity considerations. An 
approximate method of computing blade- surf ace flow properties is also derived. 

In addition, the overall pump-design problem is discussed. Methods of ob- 
taining quantities needed in the design calculations are given, and the numerical 
methods of attacking the design equations and a block diagram for a digital- 
computer program are presented. Finally, a numerical example illustrating the 
use of the method and suggestions for further use are given. 


DESIGN PROBLEM 

The general design problem can be considered from two points of view, 
namely, the engineering and the mathematical. The engineering approach is: 

Given a flow rate and a head rise, find a pump that will produce the desired 
performance with maximum efficiency. The mathematical approach is: Given the 

distribution of some flow property, such as velocity, on the boundaries (hub, 
shroud, and blade surfaces) and assumptions regarding the type of flow, find the 
impeller geometry that will result in these flow properties. 

Fundamentally, the problem is an engineering one; that is, the pump must 
produce a certain flow and a certain head rise and possibly meet other restric- 
tions, such as a maximum allowable diameter, a specified rotative speed, and so 
forth. Commonly, the engineering goal is achieved through the combined use of 
a certain amount of mathematical design, the designer's art and experience, and 
possibly development. 

The ideal approach would be to convert the engineering specifications into 
inputs for a complete mathematical design, to specify desirable flow properties 
along the impeller boundaries, and to compute the geometry of the impeller. 

Unfortunately, complete three-dimensional design methods, such as suggested 
in reference 4, are long and complex. Even when a high-speed digital computer 
is available for the computation, the setting up of the problem for the computer 
would be a long and tedious process. An exact two-dimensional method is used in 
references 5 and 6 to obtain blade-to-blade properties in a prescribed pump 
stream- tube impeller geometry without and with splitter vanes, respectively. 

This method gives an adequate picture of the flow in the leading-edge region and 
would be useful for cavitation considerations; however, the exact solution is 
quite complex and involves considerable calculations. It is important, there- 
fore, that rapid easy-to-use methods accurate enough for engineering purposes be 
developed. Two such methods of obtaining flow properties on the blade surfaces 
are given in references 1 and 7. In references 5 and 6, it was found that these 
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methods gave reasonably accurate results throughout most of the blade passage. 
Although these methods are of limited value for studying cavitation conditions, 
they are sufficient for eddy detection and for boundary-layer and loading studies 
except in the immediate vicinity of the trailing edge and in a region near the 
leading edge that extends further into the impeller as the angle of attack devi- 
ates from the design value. 

The hub-shroud design method presented herein (with optional blade-surface 
computations) is another such rapid and easy-to-use method. It enables the 
designer to proceed from a known streamline and its velocity distribution to an 
adjacent streamline and its velocity distribution. Thus, given the conditions 
along the hub, the entire hub-shroud profile is built up by proceeding from the 
hub streamline to the next streamline, and so on, until the shroud is reached. 
Conversion of the pump specifications into input for the design calculations is 
described in the DESIGN INPUTS section. 


DESIGN EQUATIONS 
Stream Tube Volume Flow 

The coordinates of an unknown streamline are related to those of an adjacent 
known streamline by a form of the continuity equation. In the derivation of this 
equation, it is assumed that surfaces of revolution obtained by rotating stream- 
lines about the impeller axis are stream surfaces and that all passages between 
blades carry the same volume flow. Then the equation for the volume flow per 
unit blade passage AQ/M through a cross section of a stream tube bounded by 
adjacent surfaces of revolution n^ and n? and by adjacent blades 0-j. and 

0^, as shown in figure 1, is obtained by 

/ r n 2 r° d(n) 

£A = / / W m r d0 dn (l) 

* J *1 \(n) 

where A is the area of a cross section of stream tube everywhere normal to the 
meridional or through-flow direction, W m is the meridional or through-flow 

velocity component and is a function of n and 0, and r is the radial dis- 
tance to the center of the element and is a function of n. (All symbols are 
defined in appendix A. ) 


Velocity Gradient 

The velocity-gradient equation (which is derived in appendix B) is used to 
obtain the velocity on a new streamline when it is known on an adjacent stream- 
line. Simplification of the exact three-dimensional equations of motion to a 
simple velocity-gradient equation is accomplished in two major steps. In the 
first step, the problem is reduced to two dimensions by considering the flow on 
a known stream surface S (fig. 2) that extends from hub to shroud throughout 
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the impeller. Mathematically, this surface can he any stream surface between 
blades. It is reasonable to assume that the mean flow follows a mean blade sur- 
face throughout the guided portion of the impeller. At the inlet, the surface 
deviates from the mean blade surface for nonzero angle of attack and, at the 
outlet, it differs because of slip. These two factors can be taken into account 
and are discussed in the section on Blade Properties. 


In general, flow properties are functions of 0 , r, and zj but on the sur- 
face, 0 Is a function of r and z. Therefore, the flow on the surface is a 
function of r and z only; that is, it is mathematically two dimensional. As 
a result, the design calculations can be set up as if the flow were in an r-z 
or meridional plane. A meridional plane is a plane that contains the Impeller 
axis, but it is not a physically meaningful plane in the flow field (except in 
the case where there is no relative circumferential velocity). An advantage of 
working in the meridional plane is that the meridional streamline picture is a 
true representation of the through flow associated with the surface S. 


In the second step, the two-dimensional equations are further simplified by 
considering flow conditions along meridional streamlines and gradients along 
normals to these streamlines (fig. 3). This consideration results in the 
velocity-gradient equation 


dW 

dn 


aW - b 


( 2 ) 



Axial distance, z 


(a) Meridional streamline. (b) Velocity components. (c) Normal to meridional stream- 

line. 

Figure 3. - Meridional plane. 


where 

cos 2 p sin 2 p cos a B cos p sin p sin a 
a — ' 1 — ■ — 

r c r r 


( 3 ) 
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and 


b = 2eo sin 3 cos a + B cos 3 


'dWg 

L dm 


+ 2a) sin a 


g ^i a) dA 
W dn + W dn 


(4) 


In equation (2), W is the fluid velocity relative to the rotating impeller and 
n is distance along a normal to the meridional streamline. All quantities in 
a and b are given as design inputs (e.g. , 03 ) or are computed from the geometry 
(e.g. , a) or from flow properties (e.g. , dW^/dn) of a known (i. e. , initially 

prescribed or previously calculated) meridional streamline. These quantities are 
discussed as they appear in the derivation in appendix B. Briefly, from the 
meridional streamline geometry shown in figure 3(a), 


dr 

tan a = — 
dz 


dm = ^dr 2 + 


da 

dm 


and from blade-shape specification, bo/br and bo/b z, 

-n bO . bo 

B = r — sm a - r — cos a 
oz or 

tan (3 = r — sin a + r _ cos a 
br b z 


Physically, (3 is the angle measured from the meridional projection to the 
streamline of the flow surface (fig. 4). The angle is positive when measured in 



Figure 4. - Projection of stream surface onto meridional plane. 
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the direction of rotation. The remaining quantities in b are obtained as 
follows: W q equals W sin |3; to , the rotative speed of the impeller, is pre- 

<3% AH; 

scribed; - — = - — where AH; is obtained from a prescribed distribution of 
’ dn An 


at the inlet and An 


is obtained from continuity; and 


dA 

dn 


AA 

An 


dHp 

as for s — . 
dn 


Hi 


With the two equations, velocity (eq. (2)) and continuity (eq. (l)), the 
hub- shroud profile can be constructed from a given meridional streamline and its 
velocity distribution and a given blade shape. The velocity on the stream 
surface throughout the impeller is of necessity obtained as part of the design 
calculation. 


Head 

The static head h along a streamline can be obtained from equation (B22): 


u W 2 , co 2 r 2 coA 

h = H i " 2i + — " T 

where Hj_ is the inlet total head that is given for each streamline. 

The total head along a streamline can be computed from equation (B28): 

„ TT o&Vq a,}, w , ovWq od 2 i* 2 o>A 

H = Hi +— - T - % +— +— " -g 

This completes the incompressible nonviscous flow calculations on the stream 
surface S. 


Blade-to-Blade Calculations 

Approximate blade-surface velocities and static heads can be obtained from 
equations (C3) to (C7). These equations are derived in appendix C with the 
assumption 6f a linear variation of static head from blade to blade. They pro- 
vide results of satisfactory accuracy except in the leading- and trailing-edge 
regions of the blades. If more exact results are required, the blade- to-blade 
analysis method of references 5 or 6 can be used. 


DESIGN INPOTS 


The equations presented in the previous section enable computation of the 
internal-flow conditions and the shroud shape of a pump impeller provided that 
the following quantities are known, that is, given or prescribed: volume flow 
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rate, Q; head rise, AH; inlet- total-head distribution, H^; inlet-pre-whirl dis- 
tribution, A; rotative speed, co; meridional contour of some streamline such as 
the hub, r = r(z); velocity along this streamline, W = W(r,z); the mean blade 
surface in terms of the curvature components, de/dr and dd/Bz; the blade thick- 
ness, t n or t q ; the number of blades, M; and any deviation of the flow surface 

from the blade surface such as slip or angle of attack. 

The flow rate Q and the head rise AH will usually be given as pump 
specifications. The rotative speed co is sometimes given. In general, co, Q, 

AH, Hj_, and A are not directly related to hub shape, velocity along the hub, 

blade shape, and number of blades, so the designer has considerable freedom in 
prescribing these quantities. At the inlet and at the outlet of the pump, how- 
ever, certain relations among these quantities and Q and AH must be satis- 
fied. 


Outlet Conditions 

At the outlet, two equations must be satisfied. The first is Euler's tur- 
bine equation: 


AH = — (rW sin 3 + cor 2 - A) av 

o 

In equation (5), q is the hydraulic efficiency defined by 


^-actual 

^■ideal 


(5) 


The hydraulic efficiency may be estimated from a consideration of losses or 
it may be assigned from experience. The rotative speed co may be chosen to be 
consistent with other free choices in equation (5) or it may be predetermined by 
other considerations; for example, the co of a direct-drive turbine that will 
drive the pump. 


If slip is to be taken into account, equation (5) becomes 


AH = 


ryo 

g 



sin p + cor) - A 

Jav 


where f g is the slip factor and is defined and discussed in the section on 
Blade Properties. 


The second equation is the equation for volume flow rate: 


Q = (W cos p) av 



(6a) 
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(6b) 


or for a radial outlet when a = 90° a.nd r s = r h 


Q = (W cos 3) av ^2*r - Mt 0 j (z s - z h ) 


av 


The values of any quantities (W, |3, etc.) needed in the design must be pre- 
scribed at the outlet to be consistent with equations (5) and (6) and with any 
other constraints imposed, for example, the existing geometry of the system in 
which the pump will be used. In equations (5) and (6) the subscript av refers 
to average values between hub and shroud. In practice, however, these equations 
are used to obtain input values on the prescribed streamline. The result of this 
procedure is that the pump may produce a somewhat different head rise or flow 
rate from that specified. If this is a significant deviation, the prescribed 
streamline and its related input can be adjusted after a trial design. This 
adjustment can be repeated as often as required. 


Inlet Conditions 

There are four possible situations that may confront the designer at the 
pump inlet: (l) The upstream flow conditions are given, for example, when an 

inducer precedes the pump; (2) the upstream geometry is given, for example, a 
predetermined inlet pipe; (3) the geometry immediately upstream of the pump is 
the responsibility of the pump designer; and (4) there are no restrictions on the 
inlet conditions, for example, it is the responsibility of someone else to make 
the inlet section fit the pump. 

For item (2), the upstream section can be analyzed (from essentially the 
same equations as those used for the pump design) and the flow conditions are 
then known. This situation is now the same as item (l). For item (3), the in- 
let section can be designed by the method of this report (with certain changes 
because there are no blades) then the flow is known as in the situation in 
item (l). Thus, for inlet conditions required as design input, there are two 
different situations, namely, items (l) and (4). In both cases, the relations 
of continuity 


AQ = (2itr - Mt 0 )An W cos 3 


(V) 


where 


r 2 - r l 
cos a-]_ 


except in the case of a radial inlet a = 90° where 


An = Zg - Z]_ 
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and conservation of angular momentum 


W sin (3 


1 


— - cor 
r 


( 8 ) 


must be satisfied at each streamline for item (l) and at some mean streamline for 
item (4). 

Note that equation (7) is written for a station just inside the impeller in- 
let rather than just upstream of the inlet and, therefore, includes blade block- 
age. Also, note that in equation (8) it is assumed that no work has yet been 
done on the fluid at this station. Not all quantities involved in equations (7) 
and (8) are actual design input at every streamline. Only the prescribed stream- 
line (usually the hub) has a W and r as input. Values of A and Hj_ (which 

does not appear in the equations) are required input at every streamline even 
when upstream conditions are not prescribed. Values of 0, tg, and M across 

the inlet are incorporated into the blade input. Angle 0 is the angle of the 
flow surface not the blade surface, so that when these angles are not equal 
(when the angle of attack is not zero) the blade shape at the inlet must be 
specified in such a way as to produce the required value of 0. 

When equations (7) and (8) are used to relate average conditions, as in 
item (4), AQ becomes Q, An becomes (r s - r^)/cos a av or z s - z^, and other 

quantities are averages from hub to shroud. Values for the prescribed streamline 
are estimated from the hub-to- shroud average values. After a design is made, the 
inlet should be reexamined at each streamline to be sure that conditions are sat- 
isfactory. In particular, the blade angle at each streamline should be such as 
to produce the expected angle of attack (usually zero). If the expected angle of 
attack is not attained, some inlet condition or the blade shape at the inlet must 
be changed. 

In cases where the upstream flow is known before the pump is designed, the 
streamlines resulting from the design calculation must match the upstream stream- 
lines at the inlet. If these two sets of streamlines do not match, some input 
condition must be changed, and a new design calculation carried out. This proc- 
ess is repeated until a satisfactory match is achieved. 


Hub Shape and Velocity Distribution 

After the inlet and outlet conditions are determined, the meridional contour 
and velocity distribution of some streamline must be chosen throughout the pump. 
This could be the hub, the shroud, or any other streamline in between. It is 
probably best to start with the hub, since streamline spacing is most sensitive 
to velocity changes near the hub and undesirable resulting shapes in succeeding 
streamlines can be more readily eliminated. Usually the hub radius r and 
velocity W are given as functions of z. Both r and W may be prescribed 
directly, or, for better control of cavitation, the static head h may be pre- 
scribed together with either W or r. These three quantities are related in 
equation (B22). 
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If static head and velocity are prescribed, the hub shape is obtained from 


r = i “\/2g(h - H i ) + W 2 + 2o>A 

If static head and hub shape are chosen, the velocity is found from 

W = V2g(Hi - h) + - 2coX 

The velocity distribution on the hub, according to air- compressor experi- 
ence, should be accelerating, if possible. At least unnecessary negative veloc- 
ity gradients should be avoided. The static-head distribution should be such 
that cavitation is avoided or controlled. No useful criterion is known for hub 
shape except perhaps ease of fabrication. 

If it is more desirable to have direct control over conditions on the 
shroud, this can be done, and the design method can be made to proceed from the 
shroud to the hub by changing the sign of An in the equations. 


Blade Properties 


When the mean flow surface is assumed to follow the mean blade surface, the 
blade shape is usually prescribed in two parts: a mean blade surface and a 

thickness distribution. In general, the mean blade surface is of the form 
9 = 0(r,z). For very high speed wheels, it is often limited to radial elements 
and is of the form 9 = 0(z). The blade surface should be prescribed in such a 
manner that the angle p can be computed conveniently at any point. One method 
would be to prescribe Ss/Sr and b9/b z, the blade- curvature components in the 
r- and z-directions, respectively. The angle p is found from the relation 

, de Se / nX 

tan p = r t— sin a + r tt— cos a (9) 

or oz 


Equation (9) results from the relation 


tan p 


d0 <d0 dr c)0 dz 

dm or dm oz dm 


which holds along a streamline where 0, r, and z are functions of m. 

For the case of radial blade elements (ho/br = 0) or axial elements 
(c)0/d z = 0), the following method may be used to determine the mean blade sur- 
face. Instead of prescribing the blade angle, either the total head rise 
H - or the relative tangential velocity along the hub is prescribed. Then 

P along the hub is found from 


P 


sin 
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where 


w 0 = 


g(H - Hi) 
cor 



cor 


For radial blade elements 


c*0 1 tan |3 

= cos a 

where p and a are known as functions of z along the huh; therefore, since 
bo/bz is independent of r, bo/bz is known everywhere as a function of z. 

For axial blade elements 


bo 1_ tan p 

3r ~ r h sin a 

and bo/br is known everywhere as a function of r, since p and a are known 
along the hub as functions of r. 

Although 6 is not needed in the design procedure, it is usually needed for 
fabrication and can be obtained from 


d 0 


_ be 

br 


dr 


* 

+ dz 
oz 


Since the flow near the outlet of the pump does not closely follow the 
blades, a closer approximation to actual flow can be obtained if a reasonable 
value of the slip factor can be estimated. The slip factor f s is defined as 
the ratio of the absolute tangential velocity of the fluid at the outlet to the 
absolute tangential velocity the fluid would have if the flow angle were equal 
to the mean blade angle, that is, if the fluid were perfectly guided by the 
blades. The effect of slip is taken into account in computing the flow angle p. 
At the outlet, p is computed from the slip factor, the mean blade angle, and 
the flow rate. In reference 7, a parabolic variation in sin p is assumed to 
hold from a point in the pump where the effect of slip begins (and where p is 
known from eq. (9)) to the outlet. Angle p is then computed in this region 
from the parabolic variation instead of equation (9). The point where slip is 
assumed to become significant can be estimated from experience or can be computed 
from an empirical formula such as that given in reference 7. 

A similar situation exists at the inlet where the actual mean stream surface 
may differ significantly from the mean blade surface in the case of nonzero angle 
of attack. The angle of attack, or incidence angle, at some streamline is 
defined as the angle between the mean flow direction at the inlet and the mean 
blade surface at the blade leading edge. Then the p of equations (7) and (8) 
is not the p of the mean blade surface. The mean blade surface at the inlet 
then has to be specified in such a way that its angle is equal to angle p plus 
the angle of attack. The effect of angle of attack extends somewhat into the 
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impeller and could "be accounted for in a manner similar to that discussed pre- 
viously for slip factor. Note that even though the mean blade angle at some mean 
streamline satisfies equations (7) and (8) at the inlet, there may still be 
nonzero angle of attack at other streamlines and, if significant, should be taken 
care of as done previously. Reference 2 discusses computation of angle of 
attack. 

The blade thickness may be prescribed in any manner that allows its compu- 
tation as a function of z and r. Normally, t r is prescribed, but the thick- 
ness in the circumferential direction tg is required for the design method and 
is given by 



The initial choice of the number of blades is best determined by experience. 
The number need not be constant, that is, splitter vanes or partial blades may 
be added at various stations throughout the pump impeller. Since the flow is 
assumed to be periodic, the number of blades (including splitter vanes) at any 
station should be an integral multiple of the number at the previous station. 

The hub shape, velocity along the hub, blade shape, and number of blades may 
be changed if the initial values result in undesirable flow conditions (see 
NUMERICAL EXAMPLE AND DISCUSSION). 


NUMERICAL PROCEDURE 

The numerical solution of equations (l) and (2) is based on the fundamental 
assumption that the distance An between adjacent streamlines is small enough 
so that properties can be assumed to be constant across An. 

In equation (2), assume that a and b are constant across a stream tube 
from np to n 2 * If equation (2) is multiplied by the integrating factor 

e /-a dn^ pp can pg integrated to give 

I W 2 = Wpe a ^ + | (1 - e a *>) (11) 

! where An = n 2 - np. 

i 

1 

Equation (ll) is used to compute W 2 when Wp is known. The parameters 
a and b are computed along the streamline at np, and An is computed from 
the next equation to be developed. 
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Equation (l) can be solved by applying the mean- value theorem. First, con- 
sider the integration with respect to 0: 



where W m is evaluated at some 0. Assume that W m is the W m of the mean 
surface S. At any value of n 


r 


je t (n) - 0 d (n)J 


2rtr - Mtg 
M 


where r and t q 


are functions of n. Now equation (l) can be written as 


AQ 



- Mtg)dn 


Assume that W m (2itr - Mtg) is constant across An. Then 


AQ = W m (2jcr - Mt 0 )(n 2 - n-jJ 


Hence 


( 12 ) 


An = 


AQ 

W m (2itr - Mtg) 


(13) 


where W m (2jtr - Mtg) is evaluated on the streamline at n^, as are the parameters 
a and h of equation (11). 

Coordinates of a new streamline are obtained by assuming that a is con- 
stant across a stream tube and by integrating equations (B13) and (B12), 
respectively ^ to give 


and 


r 2 “ r l + cos a 


Z2 = z i - An sin a 


(14) 


(15) 
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The two normal derivatives that occur in b are computed from streamline- 
input data for a linear variation with n assumed at each station: 



( H i> 2 - (H i ) 1 

(16) 

An 

\ ‘ \ 
An 

(17) 


1 


The numerical method employed to compute derivatives at discrete points 
along a streamline is a "spline curve fit’* technique taken from the SHARE Program 
for IBM 704 users. Briefly, the method requires a set of points and the slopes 
at the end points for which a cubic polynomial is defined for each interval 
(between adjacent points) such that the function represented by the set of cubic 
equations has continuous first derivatives within the given range. Since the 
end-point derivatives are usually not known, they too must be computed numeri- 
cally by a method such as Newton's end-point formulas (ref. 8). 

A block diagram for a digital- computer program is shown in figure 5. The 
diagram is explained by block number as follows: 

(l) Read pump input data: 

(a) Number of stations (points along the streamline), J 

(b) Number of streamlines, K 

(c) Rotative speed, 

(d) Flow rate, Q 

(e) Number of blades at each station, M 

(f) Prerotation for each streamline, A 

(g) Inlet total head for each streamline, Hp 

(h) Coordinates of hub for each station, z, r 

(i) Relative velocity at each station on hub, W 

(j) Tabular blade-shape data as functions of r and z, bdfbz, be/br 

(k) Tabular blade thickness as function of r and z, t n 
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For the remaining steps, the subscript j refers to stations along a streamline 
and is in the range j = 1 , 2, • . . , J. 

(2) Compute (dr/dz) ^ from spline curve fit, from equation (B8), and 

m, from 




'dr 
i dz 


11/2 


dz for 2 s j S J 


where m^ = 0 


(18) 


(3) Compute (da/dm) j from curve fit 

(4) Obtain blade-shape data, (de/dz)j and (be/dr) j, by interpolation from 

the tables (step (l)(j)) 

(5) Compute (3. from equation (9) 

(6) Compute (W g) . from equation (B18), (W m ) . from equation (B17), and 

J 3 

(dW^/dm^ from curve fit 


( 7 ) 

( 8 ) 
( 9 ) 


Obtain blade- thickness data (t n ) . by interpolation from the tables 

J 

(step (l)(k)) and compute circumferential thickness (tg)j from 
equation (10) 

Compute Bj from equation (B16) 

Compute static head hj from equation (B22) and total head Hj 
from equation (B28) 


Steps (10) and (ll) are the blade- surface calculations and are optional, since 
they are not essential steps in the design procedure. 

(10) Compute (Ah) * from equation (C3), (h,j) . from equation (C4), and (ht) . 

from equation (C5) 

(11) Compute (W(j) . from equation (C6) and (Wt) . from equation (C7) 

J J 

(12) Write out desired streamline data. In working from the hub to the 

shroud, the first streamline data will be that associated with the hub. 
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Quantities that occur in the calculation procedure of potential inter- 
est are: 


(a) Streamline coordinates 

(b) Resultant and component velocities 

(c) Meridional streamline curvature 

(d) Blade- curvature components 

(e) Angles a and p 

(f) Static head and total head 

(g) Blade- surface velocities and static head 

(h) Circumferential blade thickness 

(13) Test whether or not the solution has reached the shroud* (If the shroud 
is the starting streamline, the logic is still the same and the design 
is finished when the hub is reached; however, if some intermediate 
streamline is the initial one, the design has to proceed from this 
streamline to both the hub and the shroud, and the logic must be 
changed slightly. ) If the solution has reached the shroud, go to 
step (17); if it has not, go to step (14). 

from equation (16), 

j ' h 

from equation (17), a^ from equation (3), and bj from equation (4) 


(14) Compute (An ) * from equation (13), 


(15) Compute next streamline coordinates and velocity: r^ from equa- 

tion (14), Zj from equation (15), and Wj from equation (11) 


(1G) Return to step (2) 
(17) Scop 


NUMERICAL EXAMPLE AND DISCUSSION 

The digital- computer program as outlined in the preceding section was 
applied to a numerical example for the purpose of illustrating the use of the 
method. The following is a list of the conditions prescribed: 

(l) Number of stations, 26, approximately equally spaced with distance along 
meridional hub contour 
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Radial di 


(2) Number of streamlines , 40 

(3) Rotative speed, o>, 1571 radians/sec 

(4) Flow rate, Q, 48 cu ft /sec 

(5) Number of blades, M, 4 at each station (no splitter vanes) 

(6) Prewhirl, A, 0; inlet total head, Hj_, 1304 ft 

(7) Hub contour, r = r(z), as in figure 6(a) 

(8) Hub velocity distribution, W, shown in figure 6(b) as function of m/m 



Axial distance, ft 
(ft) Hub contour. 


Percent of distance along meridional streamline, m/mj 
(b) Hub velocity distribution. 


Figure C. - Hub input data for numerical example. 


(9) Blade-shape data, bd/bz or bd/br, as functions of z. The quantity 

bd/bz is a linear function of z from 0.000 to 0.400 as shown in the 
following table: 





Radial distance 


Axial distance 
from impeller 
inlet, 

z 


Blade- shape 
specification 


0.000 
. 050 
.150 
.250 
.350 
.390 


-10. 39374 
-9.09453 
-6. 49609 
-3.89765 
-1.29922 
-.25984 


-.12992 

0 


(10) Constant blade thickness, t n , 0.0075 ft 

Figure 7(a) shows some of the computed streamlines from 'the hub to the 
shroud. Figure 7(b) shows the velocity distributions along each of these 


Shroud 
1st trial 
2nd trial 

Streamline 
30 'l 

20 \ 1st trial 

10 J 

Hub, computed with 
shroud as input 
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(a) Computed streamlines from hub to shroud. (b) Velocity distributions along computed stream 

lines from hub to shroud . 

Figure 7. - Meridional streamlines and velocity distributions for numerical examole. 






streamlines. Figure 7 represents the first trial to obtain acceptable shroud 
shape and shroud velocity distribution with the conditions prescribed pre- 
viously. 

In order to improve the shroud velocity distribution, it was decided to 
distribute the blade curvature over a greater axial distance, in particular, to 
vary hd / Sz linearly with z to a value of z of 0. 500 instead of 0. 400. The 
new blade data is shown in the following table: 


Axi al di s t an c e 


Blade- shape 


from impeller 


specification 


inletj 

z 




3? 


0.000 
.050 
.150 
.250 
. 350 
.450 


-10. 39374 
-9. 35437 
-7.27562 
-5.19687 
-3. 11812 
-1. 03937 


0.0 


. 490 
.495 
.500 
.505 
.510 
.550 
.600 


-.20786 

-.10394 

0 


The second trial shroud shape and shroud velocity distribution are shown by 
dashed lines in figures 7(a) and (b), respectively. The number of blades, hub 
contour, and hub velocity may also be changed to improve the solution. This 
process can then be repeated until an acceptable shroud and shroud velocity dis- 
tribution are obtained. 

Although it was not done for the numerical example the blade- surface veloc- 
ities, in an actual design, also would be computed for each meridional stream- 
line, which would enable the designer to exercise a certain amount of indirect 
control over these velocities. It is usually not practical to eliminate all 
undesirable qualities of the velocities. For example, all blade-surface decel- 
eration could not be eliminated. The chief utility of the blade-surface compu- 
tations is that they can indicate severe undesirable gradients and the presence 
of an eddy on the driving surface. It is usually fairly easy to eliminate the 
eddy in later trials by raising the hub velocity, changing the blade shape, or, 
more commonly, by adding splitter vanes between blades at the station where the 
negative velocity first appears. For numerical examples of blade-surface veloc- 
ities computed by two approximate methods (including the one presented in this 
report) and their comparison with exact solutions see references 5 and 6. Ref- 
erence 6 gives results for a pump with splitter vanes. 
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Although the approximate blade- surface computations also permit indirect 
control over the blade-surface static head or pressure, they are unfortunately 
of quite limited value in studying cavitation conditions, as is pointed out in 
reference 5. If the time and effort were warranted, however, the method of 
reference 5, applied to only one stream tube, might give sufficient information 
for cavitation studies. (In ref. 5, angle of attack was not taken into account 
in the approximate calculations so that the off-design case shows worse agreement 
with the exact solution at the inlet than would be the case if angle of attack 
were included, as suggested in this report. ) 

As mentioned in the Design Inputs section, the method can also proceed from 
shroud to hub by making the sign of An in equation (13) negative and by pre- 
scribing conditions along the shroud. Adjustment of conditions on the shroud 
similar to those done to the hub can produce an acceptable hub and hub velocity 
distribution. The shroud contour and velocity distribution computed as the first 
trial (figs. 7(a) and (b)) were used as input to the program to illustrate start- 
ing on the shroud and the computed hub and hub velocity distribution, as shown 
in figures 7(a) and (b), are almost identical with the original prescribed hub 
and hub velocities. 

In choosing the number of stations for a solution, some care should be 
exercised to avoid very small spacings especially in regions of large curvature 
in order to prevent numerical difficulties that become increasingly worse as the 
solution proceeds toward the shroud. The reason is that the spacing gets smaller 
near the shroud and may become of the same order of magnitude as the error in the 
numerical procedures, or the spacing may become negative. In either case, the 
solution becomes meaningless in such a region. 

A sufixcient numuer of streainu-ines jhoold be chosen so that the resulting 
An is small enough to make the assumption of constant a and b in equa- 
tion (2) a valid one. The value of An varies with the size of the impeller 
and can usually be determined in one or two trials. It is probably not worth 
starting with less than 10 streamlines. Limited experience indicates that no 
difficulties result from having too many streamlines ; however, it is a waste of 
time to have more than necessary for the degree of accuracy desired. In the 
numerical example presented, the results of the solution from the shroud to the 
hub indicate that An was chosen small enough to make the assumption of constant 
a and b satisfactory. 


Lewis Research Center 

Rational Aeronautics and Space Administration 
Cleveland, Ohio, September 19, 1962 
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APPENDIX A 


SYMBOLS 

A area of stream tube normal to through flow 

a geometric parameter, eq, (3) 

E blade- shape psirsimsfsr, ecp (El 6) 

b geometric parameter, eq. (4) 

f s slip factor 

g acceleration due to gravity 

H absolute total head, h + V 2 /2g 

AH head rise across pump 

H re i relative total head, h + W 2 /2g 
h static head 

h* static head on surface S 

& h d - h t 

J number of stations 

K number of streamlines 

M number of blades 

m distance along meridional streamline 

n distance along normal to meridional streamline 

An normal distance between adjacent meridional streamlines 

Q volume flow rate through pump 

AQ volume flow rate through annular stream tube 

r radial distance from axis of rotation 

r c radius of curvature of meridional streamline 
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S stream- surface function and surface itself 
t time 

t blade thickness normal to mean blade surface 
tg blade thickness in circumferential direction 
V absolute fluid velocity 

W fluid velocity relative to rotating impeller 

z axial distance from impeller inlet 

a angle meridional streamline makes with impeller axis, tan" 1 dr/dz 

(3 angle between streamline of surface S and its meridional projection, 

tan" 1 r d0/dm 

tj hydraulic efficiency 

Q angular distance from radial line rotating with impeller, radians 
G t - 0 d 

A prewhirl, r ±(Vq)^ 

co rotative speed of impeller 

Subscripts: 

av average between hub and shroud 

d driving surface of blade 

h hub 

i impeller inlet 

j index for stations along meridional streamline 
m meridional component 

r radial component 

s shroud 

t trailing surface of blade 
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z axial component 
9 circumferential component 
1 known streamline 


2 


adjacent unknown streamline 


APPENDIX B 


DERIVATION OP VELOCITY - EQUATION 

The objective in this appendix is to derive a simplified relation between 
the known relative velocity on some known streamline and the unknown relative 
velocity on an adjacent streamline. Other quantities involved in the relation 
are certain properties of the known streamline and the impeller blade geometry. 

The starting point is the equations of motion governing the steady relative 
three-dimensional incompressible nonviseous flow at any point in a turbomachine 
(ref. 4): 

dW r (W 0 + mr)2 ^ Sw r W 0 Sw r ^ c>W r (W 0 + mr) 2 ^ 

dt r 53 W r 3F" + T ST + r " " s 3? (Bla) 


1 <2rV 0 dW 0 W r W 0 
r dt " dt + r 


+ 2toW r 



dW z 

dt 



(Bib) 

(Blc) 


where V 0 = W 0 + <jdt, h is static head in feet of fluid and r, Q, and z are 
cylindrical coordinates relative to the rotating impeller (see sketch). 


dr dz 



The problem can be reduced from three dimensions to two dimensions by first 
considering a relative stream surface S (fig. 2) that extends from hub to shroud 
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about midway between blades. (This is a surface of the second kind in the ter- 
minology of ref. 4. ) For steady relative flow, such a surface is a three- 
dimensional surface that rotates with the impeller and is given by a relation 
among the coordinates of the form 


S(r,0,z) ts 0 


or solving for 9 


9 = 0(r,z) 


(B2) 


This equation is used to relate the static head h 
flow field with the static head h* on the surface S. 


but on the surface S 


h a h(r,0,z) 


h* = h[r,0(r,z),z] = h*(r,z) 


of the three-dimensional 
In general, 

(B3) 


(B4) 


since 9 on the surface is not an independent variable. The relation between 
partial derivatives of the static head in the three-dimensional field with those 
on the surface S is given by 


Sh* 

Sr 

Sh* 

Sz 


Sh Sh Se 
Sr Se Sr 

Sh Sh Se 
Sz Ss Sz 


Substitution of equations (B5) in equations (Bl) yields 


(B5) 


dw r (W 0 + cur)' 


dt 


1 

r dt 


dWg 

dt 


Sh* Se i Sh 

3T + gr 3? r 3e 


W r Wg rr i Vi 


3e 


(B6a) 

(B6b) 


d ^ z Sh* S0 1 Sh 
at s SV + gr 5;733 


(B6c) 


Equations (B6) are seen to be the same as Lorenz 1 s equations for axi all y 
symmetric flow (ref. 9, p. 991) with a blade force (appendix B of ref. l). 
Instead of assuming axial symmetry, however, the equations are written for flow 
on a stream surface S. And instead of introducing the concept of blade force 
and having blade-force components in the equations, the circumferential pressure 
1 ^ h 

gradient — appears in all three equations. 
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Equation (B6b) can be used to eliminate the circumferential pressure gra- 
dient from equations (B6a) and (B6c) giving 


dw r vg ah* bell drV e\ 

dt r 53 ® dr r 3r\r dt j 


(B7a) 



dh* 

= _g ^r ■ 



(B7h) 


The remaining two equations involve known geometry, flow properties on the 
surface S, and their derivatives in the r- and z-directions only. The mathe- 
matical significance here is that the problem can be considered in the r-, z- 
or meridional plane. If attention is further limited to a streamline of the 
surface and its projection in the meridional plane, then r and z are related, 
and equations (B7) can be reduced to one equation. The projection of a stream- 
line of S onto the meridional plane is called the meridional streamline 
(fig. 4) and the velocity in the direction of this streamline is called the 
meridional or through- flow velocity ¥ m . The velocity components W r and ¥ z 

are related to each other and to W m by the angle a, which is the angle between 

the pump axis and a tangent to the meridional streamline (fig. 3). Along the 
streamline, r « r(z), and a is obtained from 


a = tan - -*- —• 
dz 


(B8) 


The velocity components and ¥, can be expressed in terms of ¥ m and 


a to yield (fig. 3) 


Differentiation yields 


W r ** W m sin a 


¥ z a W m cos a 


d ^r da, d ^m 

dt” B w bi cos a dt + dt“ sin a 

dW z TT da dW m I 

dt“ ** “ Wm Sln a dt + dt - COS a J 


(B9) 


The quantity da/dt can be eliminated from equations (B9) by making use of 
the curvature of the meridional streamline 


da 

1 s da _ dt 1 da 
r c dm ~ dm ** ¥ m dt 

dt 
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Making use of this relation and combining equations (B7) and (B9) yields 



— cos a + 
r c 



sin a - 



r 


dh* de 1 drV 0 

dr “ r c)r r dt 


(BlOa) 


^m . d ^m 

— sm a + -tv— cos a 
r~ dt 


dh* _ de 1 drV e 

dz r r dt 


(BlOb) 


The next step is to combine dh*/dr and dh*/dz and reduce equations (BIO) 
to one equation. To this end, take the directional derivative of h* in the 
direction normal to the meridional streamline, that is, in the direction 
a + 90°, and call the distance in this direction n. Then 


dh* _ dh* dr dh* dz 
dn “ dr dn dz dn 


But from figure 3 


dz 

dn 

dr 

dn 


ss cos(a + 90°) = -sin a 
= sin(a + 90°) = cos a 


and equation (Bll) can be written 

dh* dh* 


dh* 


dn- = g gr C ° S a ‘ g dz 


sm a 


(Bll) 


(B12) 

(B13) 


(B14) 


Multiply equation (BlOa) by cos a and (BlOb) by sin a and combine them to 
obtain 


dh 

dn 


Vg cos a W‘ 


■m B drV £> 
r c + r dt 


where, for convenience, B is defined as 


„ _ c)0 . bO 

B = r -r— sm a - r -r— cos a 
oz or 


(B15) 


(B16) 


and where h has been written for h* since there is no longer any need to dis- 
tinguish between them. Equation (B15) is usually called the force equation. 

The dimensions of its terms are those of acceleration or force per unit mass. 

The left side is the total force on a particle in the direction normal to the 
meridional streamline. The first term on the right is the normal component of 
the centrifugal force due to rotation about the axis of the impeller, the second 
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term is centrifugal force due to curvature of the meridional streamline, and the 
third term is the normal component of the force due to the circumferential pres- 
sure gradient. At first, it may appear that there should he no component of 
circumferential pressure gradient in the meridional plane since the circumfer- 
ential direction is normal to the meridional plane. It must be remembered, 
however, that the meridional plane is not a physical plane in which the flow 
takes place but is a computational tool that permits visualization of the through 
flow. The normal component of the circumferential pressure gradient is the com- 
ponent along a line in the flow surface S of which the normal is the meridional 
projection. In general, in moving along a line in the flow surface, 6 will 
change so that there can be a component of the circumferential pressure gradient 
along such a line. (For straight blades, i. e. , for de/d z and bdfdr both 
zero, there is no component of the circumferential pressure gradient in the 
meridional plane. ) 

The velocity components Wg and W m are related to each other and to the 

resultant velocity W through an angle p (fig. 4). This angle is often called 
the flow angle and is also the blade angle when the stream surface is parallel 
to the mean blade surface. It is the angle between a streamline of the surface 
S and its meridional projection. When the flow surface is a mean blade surface, 
P is known. 

From figure 4, it can be seen that 


W m = W cos p 


(B17) 


and 


w 0 = W Sin JB (B18) 

Making use of relations (B16), (B17), and (B18) (and certain derivatives 
and identities), equation (B15) can be written in a more useful form: 


dh cos^p 

to = 


+ 


(W sin p + cor)^ 
r 


cos a 


+ BW cos 



W sin p sin a 


2ao sin a 


(B19) 


Since an equation in W only and not h is desired, another relation 
between h and W is needed. Such a relation can be obtained from a different 
manipulation of the equations of motion, which is equivalent to writing an energy 
equation along a streamline. 
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cLi* 

Multiply equation (Bla) by ¥ r = equation (Bib) by W@ = r and 


d0 


dz 


equation (Blc) by ¥ z = — and add and combine terms to obtain 


1 d¥^ 

2 dt 


tD 2 rW r = -g ~ 


dt 


(B20) 


Integrate equation (B20) along a streamline between the pump inlet and any 
station to obtain 


| (W 2 - ¥ 2 ) - ^ (r 2 - r 2 ) = -g(h - h.) 


(B21) 


V 6 

Since W 2 = V 2 - 2YqO>t + o^r 2 and H = h + -jr- where H = total head, equa- 
tion (B21) can be written 


W 2 m 2 r 2 
h = H, -I-+ 

1 2g 2g 


cjoA 


(B22) 


where instead of r^(Vg)., A has been written, which is commonly called pre- 
whirl. The prewhirl is a function of the upstream conditions and is usually 
specified for the pump designer. 

Take the directional derivative of h in the direction of n, that is, the 
direction normal to the meridional streamline to obtain 


^i ¥ d¥ m 2 r dr dA 

dn - dn g dn g dn g dn 


dh 


(B23) 


Equation (B23) can be used to eliminate dh/dn from equation (B19) with 
the result that 


d¥ _ „ 
dn \ r 


cosZp sin^p cos a B cos 8 sin g sin 


a 


- 2u> sin 3 cos a 


- B cos p 


fdWg \ <3Hi 

+ 2cu sin a + 


idm 


¥ dn 


(D dA 
¥ dn 


(B24) 
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For convenience, define the quantities a and b so that 


a = 


COS 2 p 


sirrp cos a 
r 


B cos p sin p sin a 
r 


(B25) 


and 


fm 


i ~ J 9 

b » 2co sin p cos a + B cos pi- h 2co sin a 


dm 


< b26 > 


Put these into equation (B24) to obtain the velocity equation 

aw 


an 


= aW - b 


(B27) 


This equation together with a suitable form of the continuity equation 
enables the aesigner to compute the coorainates of any streamline of the surface 
S ana its velocity attribution proviaea that the stream surface is known or 
prescribe! an! some streamline of the surface together with its velocity at- 
tribution is known or prescribe! (boun!ary con!ition). Furthermore, the static 
hea! on the surface can be compute! from equation (B22). The total hea! H can 

V 2 

be foun! by putting h = H - -5— in equation (B22) which gives, after some manip- 

2 S 

ulation. 


cnrV, 


H = % + 


g g 1 


co 2 r 2 o>A 

g + ~ i s’ 


(B28) 
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APPENDIX C 


EQUATION FOR BLADE SURFACE VELOCITY AND HEAD RISE 


With further simplifying assumptions, the approximate velocity and head on 
the blade surfaces can be computed. The static-head derivative in the circum- 
ferential direction, that is, from blade to blade, is given by equation (Bib); 


1 ^0 

r dt 


dW 0 W r W 0 Sw 0 W 0 dw q c»W 0 

dt~ + ~ r“ + 2cSX r = — ^e~ + W z ^r 


4* 


w r w 0 


+ 2u>W r ss 


g dh 
r 


Making use of relations among velocity components and angles, the previous 
equation becomes 


W cos pi 


(dWg 

dm 


4* — sin J3 sin a + 2oo sin a I - 


-g dh 
r 30 


(Cl) 


If a linear variation in static head from blade to blade in the circumfer- 
ential direction (i. e. , along a path of constant z and r) is assumed. 


1 dh _ -Ah 
r <50 “ r A0 


(C2) 


where 


and 


-Ah = h t - h d 


r A9 = r(0 t - 0 d ) = 


2rtr - Mt 0 


M 


where t 0 is the blade thickness in the circumferential direction. 
Combine equations (Cl) and (C2) and solve for Ah to obtain 


d /2jtr - Mt 0 i 
^ - - 1 — — 


W cos p| 


fdWi 


,dm 


0 W . n . n 

+ — sin p sin a + 2 cd sin a 


(C3) 


If the static head h (eq. (B22)) evaluated on the surface S is assumed 
to be the average between blades. 


h = 


h d + h t 


, Ah . Ah 
h d - — = h t + — 
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which gives 


. , Ah 

hd = h + — 


and 


, . Ah 

h t = h - — 


(C4) 


(C5) 


The relative velocity on the blade surface can be computed by making use of 
the relative total head, H red : 


H t 


= h + 


W 


rel “ “ ' 2g 
Put this relation in equation (B22) to obtain 


®rel - Hp + 


<r> 2 r 2 

2g 


oik 


If H d and A are not functions of 0, at constant r, H red is constant 


and 




h t 


+ 


2g 


so that 

W d = V w2 - g Ah (C6) 

and 

W t = V y2 + S Ah (C7) 

In equation (C6), if W 2 - g Ah is negative, W d is imaginary; however, 

% is interpreted as being negative, which indicates an eddy on the driving 

surface of the blade. In reference 6, results obtained by this method were com- 
pared with an exact blade- to-blade solution of a case with an eddy. The present 
method (which interprets the imaginary W d as negative) did indicate the eddy, 

but it somewhat exaggerated the magnitude of the negative velocities. 
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